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Abstract 

Canonically, the quantum electrodynamic radiative corrections in bound systems 
have been evaluated in photon energy regularization, i.e. using a noncovariant over- 
lapping parameter that separates the high-energy relativistic scales of the virtual 
quanta from the nonrelativistic domain. Here, we calculate the higher-order correc- 
tions to the one-photon self-energy calculation with three different overlapping pa- 
rameters (photon energy, photon mass and dimensional regularization) and demon- 
strate the reparameterization invariance of nonrelativistic quantum electrodynamics 
(NRQED) using this particular example. We also present new techniques for the cal- 
culation of the low-energy part of this correction, which lead to results for the Lamb 
shift of highly excited states that are important for high-precision spectroscopy. 
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Introduction. — In 1986, Caswell and Lepage introduced the idea of nonrela- 
tivistic quantum electrodynamics (NRQED) to describe bound states [1]. The 
underlying notion is to reformulate the predictions of full relativistic quan- 
tum electrodynamics in terms of operators acting on nonrelativistic spinors, 
so that the higher-order corrections can be expressed in terms of a series of 
iterated operators of lower order, and additional operators which come in at 
every given order in the so-called Za expansion, where Z is the nuclear charge 
number, and a is the fine-structure constant. 

The development of NRQED seems to have been motivated by the well-known 
fact that a "usual" quantum field theory based on S'-matrix elements evaluated 
on free states and free propagators is not enough to describe bound states. 
Two energy scales are present in the problem, and these have to be separated 
by a so-called overlapping parameter (see § 123 of Ref. [2]). The two energy 
domains are the high-energy relativistic scale of the virtual quanta and the 
nonrelativistic domain of bound-state momenta and energies. 

Note that the regularization and renormalization of amplitudes in the ultra- 
violet (UV) at some mass scale Auv has got nothing to do with the scale- 
separation, or overlapping, parameter. For the overlapping parameter, one 
can has use either a photon energy e, or a photon mass /i, or one can work 
in A — 2eD space-time dimensions (we emphasize that e ^ Ed, the index D is 
supplied in this work in order to facilitate the distinction of the two regulariza- 
tions). Feynman [5], as well as French and Weisskopf [4J did their calculations 
in photon energy regularization. It turned out to be nontrivial to convert the 
high-energy part from a photon mass yU to a photon energy e (see the rather 
well-known footnote 13 on p. 777 of Ref. [3j), while the photon energy regu- 
larization is the most natural cutoff for the low-energy part. 

All three regularizations have been used in the literature for the treatment of 
different bound-state QED problems. E.g., Nio and Kinoshita [5] used photon 
mass regularization for their calculation of the higher-order binding correc- 
tions to the muonium hyperfine structure, while Pachucki [6] used photon 
energy regularization for the same problem. Dimensional regularization has 
been described for the lowest-order Lamb shift by Pineda and Soto in 1998 [7|. 
Higher-order binding corrections to the Lamb shift have been evaluated in di- 
mensional regularization in [8]. 

This Letter has a twofold purpose, (i) Here, for the first time to the best of 
our knowledge, a calculation of a nontrivial QED correction is presented in 
all three common (re-)parameterizations of NRQED (photon energy, photon 
mass and dimensional regularization). Namely, we consider the higher-order 
binding corrections to the one-loop self-energy in hydrogenlike systems. We 
thereby verify the reparameterization invariance to the full extent, for all three 
common regularization methods, while working on the same problem employ- 
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ing three different metliods. (ii) As a byproduct, wc develop metliods to do 
tlie calculation of relativistic Bethc logarithms for arbitrary Rydbcrg states of 
hydrogenlike systems, and we evaluate these corrections for states with prin- 
cipal quantum numbers as high as n = 12, where the excessive number of 
terms has been prohibiting both numerical as well as analytic approaches in 
the past. 

Reparameterization invariance (general remarks). — Our final goal is to ex- 
press the self-energy shift AE{nLj) of a general hydrogenic state with orbital 
angular momentum quantum number L >2, total angular momentum j and 
principal quantum number n, 

AE{nL,) = ^ F(nL,) , (1) 

in terms of the reference state quantum numbers (we use natural units with 
h — c — eo — 1, and we choose the energy scale so that the electron mass 
m — 1). The scaled self-energy function F{nLj) has the following semi-analytic 
expansion (it is not analytic because of the presence of logarithms) , 

F{nLj) = ^40 + {Zaf {a^i ln[(Za)-2] + A^o} , (2) 

where the first index of the A coefficients counts the number of factors Za, 
whereas the second counts the power of the logarithm ln[(ZQ;)~^]. 

The reparameterization invariance of NRQED implies that the F function 
should be expressible as the sum of a regularized high-energy part Fh and a 
regularized low-energy part Fl, where Fh and Fl can be formulated in photon 
energy, photon mass or in dimensional regularization, as follows, 

F = Fnie) + Fi(e) = Fniii) + F^i^) = FnisD) + Fl(£d) • (3) 



High-energy part. — In the treatment of the one-loop self-energy, we start with 
the high-energy part, which corresponds to photon energies of the order of 
the electron mass, and electron momenta of the order of the atomic momenta 
Za, where Z is the nuclear charge, and a is the fine-structure constant. We 
identify all operators that contribute at the order ai^Za)^, and evaluate these 
for a general state in a hydrogenlike system with orbital angular momentum 
quantum number L > 2, in photon energy, photon mass, and also in dimen- 
sional regularization. We find that the final expressions simplifies considerably 
for these states, and indeed very compact final results can be indicated. 

The different contributions to the high-energy part, for states with novanish- 
ing orbital angular momentum, can be described as follows, in terms of the 
electron Dirac form factor Fi and the electron magnetic form factor F2. Here, 
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we give only an indication of these corrections, detailed formulas correspond- 
ing to the terms mentioned below can be found in Ref. [8]. First, we have an 
F[{0) correction evaluated on the relativistic wave function, where the latter 
has to be expanded up to the relative order of {Za^. This correction can be 
rewritten as the sum of several effective operators acting on the nonrelativistic 
wave function. Then, we have an -^"(O) correction evaluated on the nonrel- 
ativistic wave function. From the magnetic form factor, we have an F2{0) 
correction evaluated on the relativistic wave function, and an -^2(0) correction 
on its nonrelativistic counterpart. The form factors are known in both photon 
mass piTU] as well as dimensional PDIT^ regularization. For dimensional reg- 
ularization, all the relevant formulas are explicitly given in [8], and the terms 
corresponding to the above mentioned form factors are indicated. In order to 
go over to photon energy regularization, one has to convert the photon mass 
to a noncovariant cutoff. For the lowest-order form factor slope contributing 
to the leading a [Zaf correction to the Lamb shift (in units of the electron 
mass), this is described in the textbook of Itzykson and Zuber |I3]. For a 
general hydrogenic state, we use a different ansatz, namely a generalization 
of the approach described previously for P and D states in Refs. p^fTS] , in 
order to express the high-energy part as a function of e for a general state of 
the hydrogen atom. 



There is a further two-vertex operator which is given by the diagrams in 
Fig. 5 of Ref. It corresponds to the following Hamiltonians in the three 
regularizations, 
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The above formulas, however, are of little use for a comparison to experi- 
ments unless complemented by their evaluation on a general hydrogenic state 
in terms of actual quantum numbers. A general result for the high-energy part 
in dimensional regularization, valid for all states with nonvanishing angular 
momentum and for the weighted difference of nS states (where n is the prin- 
cipal quantum number), has been given in Ref. [S]. We here refer to Eq. (3.35) 
ibid., with partial results given in Eqs. (3.32) and (3.34) ibid., and the latter 
term corresponds to our H{en). This result is expressed in terms of matrix 
elements to be evaluated on the reference state, which is manifestly taken as a 
nonrelativistic Schrodinger eigenstate. These matrix elements, as given in 
constitute rather complicated expressions and are not evaluated in terms of 
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quantum numbers. In photon energy regularization, the general form of the 
result for the high-energy part has been indicated in Eq. (8) of Ref. [15], but 
the quantities K, and C in that equation were given in general form only for 
selected submanifolds of states. 



In this Letter, we are in the position to note that the final results for the high- 
energy part, in all three regularizations, can be expressed in a very compact 
form for all states with orbital angular momentum L > 2, 
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The Aqi coefficient is defined in Eq. ([2]) and can be given as (L > 2) 
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where |(^) is the Schrodinger eigenstate. Note that Aqi is independent of j for 
L >2. The matrix element H is derived from the magnetic form factor correc- 
tion to the Lamb shift and can be expressed either as a sum of various effective 
operators acting on the nonrelativistic hydrogenic wave function, or as a sin- 
gle operator acting on the full relativistic Dirac wave function, appropriately 
expanded in powers of Za. The latter approach leads to the most compact 
expression, and the resulting matrix element can be related to the integral 
denoted as C~^„^ on p. 4483 of [16] and evaluated using generalized virial 
relations for the Dirac equation. Indeed, the result reads, expanded in sub- 
leading order in the Za-expansion [E = —W is the electric field generated 
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by the atomic nucleus with V = 
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where ip is the relativistic Dirac wave function, and ip^ , for clarity, is its adjoint 
(row vector in spinor space, complex conjugated), which is different from the 
Dirac adjoint %[) = ip^'~^^ . The Dirac quantum number is k = 2{L — j){j + |). 
We here define S2 to be the coefficient of the {Za)"^ term (this convention will 
be useful later). This completes our treatment of high-energy photons. 

Low-energy part. — In a certain sense, the photon energy regularization consti- 
tutes the most natural procedure for low-energy photons. One simply expands 
the transition current via a Foldy-Wouthuysen transformation [1^, and then 
one applies time-independent perturbation theory from the low-energy terms 
in the resulting NRQED Hamiltonian. One then integrates the photon energy 
to some upper cutoff e (in [T7j, it is explained why the expansion first in a, 
then in e is actually an expansion for large e). 

We now describe briefly how to convert the result obtained in photon energy 
regularization to photon mass regularization. For the leading-order term of 
order a{ZaY, the by now famous substitution [SPP^] reads ln(/i) — > ln(2e) + | 
while for the higher-order terms, one has to be very careful in distinguishing 

k = \k\ from u = \JhP' + 11^. The so-called quadrupole term obtained by 

— * 

expanding the exponential exp(ifc • r) in the nonrelativistic transition current 
exp(iA;-r) is very sensitive to the changes in the matching of and e because 
the power of the photon momentum k is different from the nonrelativistic 
dipole term. The additional terms can, however, be written in closed analytic 
form. 



Finally, the full evaluation of the low-energy part in dimensional regularization 
is described in detail in Ref . |Sj , and we are now in the position to indicate the 
results as follows. We denote the (nonrelativistic) Bethe logarithm by In /cq and 
the relativistic Bethe logarithm by £, following the conventions of Ref. p3][8] . 
Both of these quantities are of course state dependent, and they can both be 
evaluated only numerically. In the three different regularizations, the results 
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read (for states with angular momenta L > 2) 
4 



Fi(e) = --\nko + iZay{A, 
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Table 1 

Explicit high- and low-energy parts for the 8-D3/2 state. The Fh is the contribution 
to the self-energy correction from the high-energy part, the Fl is the low-energy 
part, and the three regularizations are: e denotes the photon energy, fi denotes the 
photon mass, and in dimensional regularization, we work in 4 — 2e£) space-time 
dimensions. 
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Adding the high- and low-energy parts. — It is easy to see that when adding the 
high- and low-energy parts from Eqs. ([5]) and ([8]), not only the regularization 
parameters cancel, but also, a reparameterization-invariant result is obtained. 



2/t(2L + 1) 
+ {Zaf I A 



- In ko 
3 " 



61 



ln(i(Za)-2) + 



17 
15 



+ ^2 + jC 



(9) 



The reparameterization invariance of NRQED is thus verified in a nontrivial 
calculation beyond leading order, in all three common regularization methods. 
A concrete numerical example is given in Table [H where the explicit numer- 
ical coefficients are written out for the 8D3/2 state (this hydrogenic level is 
spectroscopically important [TSj). 



Having obtained compact expressions, the question can be asked whether it 
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is possible to evaluate, beyond leading order, the relativistic Bethe logarithms 
L for highly excited states of hydrogenlike atoms, in approximately the same 
way as for the nonrelativistic counterparts (the "usual" Bethe logarithms), 
for which a systematic investigation has been started in Ref. [12] in relation 
to excited states. In order to appreciate the difficulties associated with the 
problem, one should recall that the relativistic Bethe logarithms represent a 
comparatively much more demanding calculation as compared to their non- 
relativistic counterparts, and the first such evaluation was not done until 1993 
(see Ref. [20]), i.e. 46 years after the evaluation of the nonrelativistic counter- 
part [21] • 

Analytic and semi-analytic calculations, where all expressions are kept in full 
analytic form before the final photon energy integration, are prohibitively dif- 
ficult for states with higher principal quantum numbers, as already described 
in a number of previous works on the subject of interest. It is doubtful if the 
analytic approach to the evaluation of matrix elements with the hydrogenic 
propagator, which is commonly based on a Sturmian decomposition [2^23f24j . 
can ever be generalized beyond principal quantum number n = 8, where on 
the order of 10^ terms are encountered in intermediate steps [15]. Calculations 
for the relativistic corrections to higher excited states seem to be possible only 
via completely numerical (lattice) methods. 

Table 2 

Relativistic Bethe Logarithms L and 
Ago coefficients for highly excited D 



states. 
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0.025 043 91(5) 


0.022 564 66(5) 


10 


0.025 185 92(5) 


0.022 669 65(5) 
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10 


0.008 413 79(5) 


0.034 832 71(5) 


11 


0.008 681 09(5) 


0.034 876 38(5) 


12 


0.008 909 60(5) 


0.034 896 67(5) 



Here, a numerical approach inspired by a discretized space as used by Salomon- 
son and Oester [25] is used, and up to eleven-point discretized representations 
are used in order to represent differential operators on the lattice whose coor- 
dinates are chosen to represent very accurately the origin in coordinate space. 
Values for the relativistic Bethe logarithms C and for the A^q coefficients of 
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highly excited D states are given in Table [21 where we note that the I2D3/2 
and I2D5/2 states are of particular experimental interest [26]. 

Conclusions. — In summary, we have completed two goals in this Letter, (i) 
The reparameterization invariance of NRQED has been verified through rel- 
ative order (Za)^ for a rather fundamentally important QED correction to 
the spectrum of hydrogenlike atoms: namely, the one-photon self-energy for 
excited states in a hydrogenlike system with orbital angular momentum quan- 
tum number L > 2. It has been verified that the photon energy, the photon 
mass and the dimensional regularizations give the same results for the energy 
shift [see Eqs. dS]), dHD and (^]. Because the higher-order binding corrections 
to the Lamb shift involve a multitude of terms, this fact is rather nontrivial 
and is displayed in a particularly clear manner in the compact expressions 
for the self-energy effects obtained here, (ii) Numerical techniques for the 
calculation of the relativistic Bethe logarithm C have been developed which 
circumvent problems associated to the growth of the number of terms in inter- 
mediate steps with the principal quantum number; these problems otherwise 
prohibit analytic and semi-analytic evaluations for highly excited states. With 
the methods described here, calculations become possible for Rydberg states 
of the hydrogen atom, and these are important for ultra-high-precision spec- 
troscopy [T5]|2^ . 

The two above mentioned aspects are important for two rather diverse topics: 

(i) for a fundamental reassurance regarding the internal consistency of NRQED 
and the consistency of overlapping parameters used in field theories in general 

(ii) for obtaining improved theoretical predictions for transition frequencies in 
hydrogenlike atoms. 
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